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Abstract. - Iron, cobalt and nickel nanoparticles, grown in the gas phase, are known to arrange in 
chains and bracelet-like rings due to the long-range dipolar interaction between the ferromagnetic 
(or super-paramagnetic) particles. We investigate the dynamics and thermodynamics of such 
magnetic dipolar nanoparticles for low densities using molecular dynamics simulations and analyze 
the influence of temperature and external magnetic fields on two and three dimensional systems. 
The obtained phase diagrams can be understood by using simple energetic arguments. 



G 

O 

o 



> 
o 

0\ 

m 

o 
o 



X 



The formation, structure and properties of nanoparti- 
cles grown in the gas or liquid phase belongs to an ac- 
tive field of basic research, which is of interest to future 
applications in different areas, from computer technology 
to catalytic and biomedical applications (e.g., see [1]). 
Magnetic nanoparticles show a variety of unusual prop- 
erties when compared to non-magnetic particles (for an 
overview, see [2]), which can often be attributed to the 
anisotropic and long-range dipolar interaction between 
the particles. For example, ferromagnetic and super- 
paramagnetic particles with low mobility can be influenced 
by an external magnetic field, leading to various arrange- 
ments and lattice structures both in experiments [3,4] and 
in computer simulations [5,6]. On the other hand, in 
recent experiments on magnetic nanoparticles with high 
mobility the formation of chains, rings and network-like 
structures has been observed [7-12]. However, the influ- 
ence of particle size, magnetic field, and temperature on 
the chain formation is not well understood and is the topic 
of this work. 

We consider a magnetic gas of N identical, homoge- 
neously magnetized spherical particles with diameter a, 
mass m, and magnetic moment pi — fip^ located at po- 
sition fi. The particles are either enclosed in a three di- 
mensional (3d) cube or restricted to a two dimensional 
(2c?) substrate, both with linear size L = Lja. We use 
fixed boundary conditions, i.e. the particles are reflected 
by the boundaries upon contact. The interaction poten- 
tial of two particles i and j at distance = fj — f is 
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Uf, , where Uf, is the long-range anisotropic 



dipolar interaction, 
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Pi ■ % ~ 3 [Pi ■ fij) (fij (1) 



with = f^ j rij , and Uh is a sufficiently rigid isotropic 
hard sphere interaction. Within this work we chose 



167T<7 3 



(2) 



which leads for TV — 2 particles to a "head to tail" ground 
state configuration with energy — Eq, where 



En = 



27TO- 3 



(3) 



The total potential energy of N particles in an external 
magnetic field B is thus given by 
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B ■ pi. 
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In the molecular dynamics simulation, we numerically 
solve the equations of motion for the positions, f{t), and 
the orientations of the magnetic moments, fii(t), of the 
particles under the influence of the force F{ = F(f, t) and 
the local magnetic field Bi = B(f,t) given by 



Fi = 



dU_ 

df ' 



B, = 



dU 
dpi' 



(5) 
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Fig. 1: (Color online): Kinetic energy e km , potential energy tt, 
and total energy e tot of a 3d system with N = 9 particles as 
function of temperature T for different system sizes L, as ob- 
tained from the molecular dynamics simulation. The transition 
temperature T c is shown as dotted line. 

As the magnetic moment fa is assumed to be pinned 
within the particle by intrinsic anisotropics, the local mag- 
netic field Bi will provide a torque Mj = //, x Bi to the 
particle. The resulting equations of motion have the form 

Pi = mr. t = Fi, (6) 
Li = Iuji = Mi = jli X Bi, (7) 
A*, = WjX fa, (8) 

were pi is the momentum, / = ma 2 /10 the moment of 
inertia, and LUi the angular velocity of the zth particle. 

The equations of motion are integrated using the veloc- 
ity Verlet algorithm [13]. The molecular dynamics simula- 
tions can be performed within the micro-canonical ensem- 
ble with constant total energy, E tot = E kia +U = const., 
and also within the canonical ensemble with constant ther- 
mal energy, 

E theim = ^Nk B T = (E kin ) = const., (9) 

where each particle has / — 6 degrees of freedom. In this 
case the kinetic energy E km is adjusted using the Andersen 
thermostat [14]. 

In the following, we present results of canonical simu- 
lations at fixed volume V, particle number N, tempera- 
ture T, and magnetic field B. We restrict N to rather 
small values below 50 particles in order to keep the sta- 
tistical errors small. As we want to compare our results 
to experiments done in the gas phase [9, 11], we consider 
systems with very low particle densities and focus on the 
case L — > oo at fixed N. Figure 1 shows the reduced 
kinetic energy e = (E km )/EoN, the reduced poten- 
tial energy u = (U)/E N and the reduced total energy 



Fig. 2: (Color online): Finite-size scaling plot of the potential 
energy u and specific heat cv (inset) as function of the scaling 
variable x. 

e tot _ (E^^/EqN versus reduced temperature for a sys- 
tem consisting of N = 9 particles, with different linear 
size L = L/a in three dimensions. N = 9 is just one 
particular choice, N > 9 yields similar results but larger 
statistical errors. It can be seen that the potential energy 
u develops a jump with growing system size L, which re- 
sults in a jump in the total energy e tot , as the kinetic part 
is simply proportional to temperature in the canonical en- 
semble (see eq. (9)). The L-dependence of u stems from 
the fact that in a finite simulation volume particles which 
leave the bound state, chain or ring, are reflected back to- 
wards the bound structure by the boundaries. Thus the 
bound state is more stable in small volumes than in large 
or even infinite volume. This effect broadens the width of 
the transition and shifts the effective transition tempera- 
ture T C (L) to higher temperatures for smaller L. 

In order to determine the nature of the transition as well 
as the critical temperature T c , a finite-size scaling plot of 
the potential energy u(T) is shown in fig. 2, together with 
the specific heat at constant volume and constant number 
of particles, cy = du/dT\y,N- Using the scaling variable 

* = GH £+ r (10) 

a data collapse both in u(T) and in cy(T) can be achieved, 
leading to an estimation for the critical temperature T c , 

k B T c = 6.0(6) x 10- 2 £ - (11) 

The constant a describes corrections to scaling and has the 
value a = 25(3). The fact that no rescaling is necessary in 
u shows that the transition is of first order. 

In the following discussion we neglect quantum mechan- 
ical effects which lead to, e.g., size-dependent magnetic 
moments in Fe particles with less than 1000 atoms [15]. 
In our model the ground state energy Eq depends on the 
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Fig. 3: (Color online): Phase diagram of magnetic nanoparti- 
cles as function of particle size a and temperature T. Shown 
are the melting temperature (dotted black), the Curie temper- 
ature (dashed blue), and the chain building temperature, both 
with constant moment \i (long dashed pink) as well as with 
temperature dependent moment /x(T) (thick red), see text. 



size a and magnetic moment fi of the particle via eq. (3), 
leading to a simple dependence of the critical temperature 
T c on the considered type of nanoparticles. For instance, 
Fe nanoparticles with a diameter of a f=a 6 nm (10 4 atoms) 
have a saturation magnetization of /i ~ 2.2 x 10 4 /ib, where 
/ib is the Bohr magneton. The critical temperature for the 
Fe nanoparticles follows from eq. (11) to T c ps 140 K. A 
resulting phase diagram for magnetic nanoparticles is de- 
picted in fig. 3, where T c (a) is shown as a dashed line. 
As to lowest order, T c {a) oc er 3 , it has to be compared 
to other characteristic temperatures like the Curie tem- 
perature Tcuric and structural transition temperatures, 
e.g., the melting temperature TmcU, in order to deter- 
mine the validity of the theory. Hence these two tempera- 
tures, which are size dependent as well, are also depicted 
schematically for typical magnetic systems. The finite-size 
form of TMeit can be approximated to lowest order using 
energetic arguments to give T Me i t (a)/T^ elt = 1 - a M cit/cr 
with material specific constant <7Moit [16]. In the case 
of the magnetic transition, standard finite-size scaling 
theory gives, again to lowest order, a temperature shift 
of the form T Cmie {a) /Tg> mie = 1 - (CT Curic / ( x) 1/l ', with 
the exponent of magnetic correlations v and constant 
ccurie [17]. As the magnetization /i vanishes at Tc ur i e 
as /x(T) = /i(0)(l — T/Tcurio)' 3 with critical exponent /3, 
and as E is proportional to /x 2 (eq. (3)), we can calculate 
a corrected T* c (a) with temperature dependent magnetic 
moments as solution of the self-consistency equation 

3™)". (12) 
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Fig. 4: (Color online): Magnetization m(B) of N — 25 particles 
confined to a 2d surface as a function of the magnetic field B for 
different temperatures T. The arrows indicate the direction of 
the hysteresis loop. The corresponding mean reduced potential 
energy u(B) is shown in the inset. 

shown in fig. 3 as thick red line [18]. For small par- 
ticles T c sc (er) w T c ((t), as the rhs. of eq. (12) is ap- 
proximately unity. On the other hand, for big particles 
T c sc (ct) « Tcuricfff), as now the lhs. of eq. (12) is vanish- 
ing. As a result, in the case of Fe nanoparticles, a phase 
transition from ferromagnetic chains/rings to a ferromag- 
netic gas can occur for particle sizes below approximately 
20 nm, while for larger particles the transition goes from 
ferromagnetic chains/rings to a paramagnetic gas. 

In the following we consider two and three dimensional 
systems in an external magnetic field B, which in the 2d 
case is oriented perpendicular to the plane. The field 
B = \B\ allows us to analyze the stability of the struc- 
tures which are observed in the simulations - namely a 
closed ring, a chain and a lattice-like structure in the 2d 
case. Typical hysteresis loops of a 2d system for different 
temperatures are shown in fig. 4. At low temperatures, 
T = 0.06 T c , the system shows a broad hysteresis loop. 
Starting from the initial ring structure with the zero field 
value, ttl — 0, the magnetization grows linearly while the 
moments of the particles start to turn into the direction of 
the increasing field. The following jump to the saturation 
value is related to the dissolving of the ring, because af- 
terwards the particles are unbonded, forming a lattice-like 
structure on the 2d surface. Then, the moments of the 
particles point into the direction of the field, which gives 
rise to the saturation magnetization of the system. With 
decreasing field, the magnetization remains saturated and 
the particles stay in the lattice-like structure until this 
structure becomes unstable due to thermal fluctuations. 
The jump down to the linear part of the magnetization 
follows from the structural rearrangement of the parti- 
cles back to the ring structure. With increasing temper- 
ature, the width of the hysteresis loop becomes smaller 
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and vanishes at T ss 0.69 T c . In the temperature range 
0-69 T c <T <T C the potential energy u(B) shows a cusp 
at B c (see inset of fig. 4), which vanishes at the critical 
temperature T c . Above T c the system is paramagnetic. 

In the following, a simple model is presented, which, 
based on analytical ground state arguments like the poten- 
tial energy, is able to explain qualitatively the numerical 
results in both 2d and 3d systems. For a large number of 
particles, N — > oo, both the ring and the chain structures 
can be treated equally regarding their potential energy u. 
Due to the long-range nature of the dipolar interaction, 
the corresponding result for the reduced potential energy 
in the ground state is u = — C(3), where £(s) denotes the 
Riemann Zeta function. In both 2d and 3d systems at 
T <T C and B — the particles form a ring, which is the 
most favorable state for N > 4 [19]. In the presence of a 
magnetic field oriented perpendicular to the system in 2d, 
the equilibrium state is changed: In the 3d case, the ring 
plane turns perpendicular to the field, similar to the be- 
havior of an antiferromagnetic spin systems in an external 
magnetic field. Furthermore, in both cases the magnetic 
moment /I, of each particle turns into the field direction by 
an angle a(B), leading to a potential energy per particle, 



u(a) = -C(3) 



^cos 2 (a) - i + V6^- sin( 



a 



(13) 



with B c = ^/6(^(3)Eo/ \x. The equilibrium state is charac- 
terized by the angle 



a(B) 




(14) 



which follows from eq. (13) by requiring du/da = 0. Using 
this dependency, we can express the potential energy u as 
a function of the magnetic field to get 



u(B) = -C(3) 



1+ — 



(15) 



This curve is shown as thin black line in the inset of fig. 4. 

The ring structure remains stable as long as B < B c 
and the dipolar interaction is attractive. This changes 
at B = B Cl where the angle a approaches the critical 
value a(B c ) = arcsin( y / 2/3) « 54.74°. Here, the dipolar 
interaction becomes repulsive, leading to a break up of the 
ring due to the mobility of the particles (note that B c is 
independent of the size of the particles; for saturated Fe 
spheres we find B c w 0.47 T). In the 3d case, the particles 
rearrange in chains aligned parallel to the field, which is 
not possible in the 2d the particles are confined to 

the surface. Instead, they rearrange in a hexagonal lattice- 
like structure on the surface, which is stabilized due to the 
repulsive character of the dipolar interaction. Now, both 
in three and in two dimensions, u is mainly dominated 
by the field energy being proportional to — B (see inset in 
fig. 4). 
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Fig. 5: (Color online) Phase diagram of the 2d system. The 
arrows indicate the different structures of nanoparticles when 
crossing the phase boundaries. Stars mark the transition from 
the bonded to the unbonded phase, triangles the reverse tran- 
sitions, as obtained from the molecular dynamics simulation. 
For clarity, the simulation data is only shown in the lower half 
plane. The dashed line marks the transition to the lattice- like 
structure with decreasing temperature. 

The reduced magnetization m parallel to the magnetic 
field can be calculated for ring structures as a function of 
the field B, which yields 



m(a) = 
m(B) = 



M(a) 
2 B 



sin(a), 



(16) 
(17) 



where M s is the saturation magnetization (thin black line 
in fig. 4). 

These ground state arguments can be expanded to de- 
scribe the influence of temperature on the system, by 
specifying the transitions between the different states by 
purely energetic arguments. At finite temperatures and 
non-vanishing magnetic fields, the ring can be broken up 
by overcoming an energy barrier Au by thermal activation. 
This energy barrier separates the ring and the lattice-like 
structure in 2c? and the ring and the chain structure in 
3d respectively. Basically the barrier height is given by 
Au(B) = u(B c ) — u(B), which is the energy difference 
between the unstable ring and ring in equilibrium in the 
presence of the field B, (eq. (13)). If the mean kinetic en- 
ergy equals this barrier height, the ring becomes unstable 
and dissolves. This and the equivalent consideration for 
the transition back to the ring structure leads to 



Bl(T) = B c 




(18) 



If the field exceeds the critical value Bl(T), the ring dis- 
solves into the lattice-like structure due to the kinetic fluc- 
tuations and the influence of the field. For fields lower 
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Fig. 6: (Color online) Phase diagram of the 3d system. The 
arrows indicate the different structures of nanoparticles when 
crossing the phase boundaries. Stars mark the transition from 
the ring to the chain phase, triangles the reverse transitions, as 
obtained from the molecular dynamics simulation. For clarity, 
the simulation data is only shown in the lower half plane. 

than 

*i(T) = *(l + V ^£ c , (19) 

the dipolar particles flip back to the ring structure. Fig- 
ure 5 shows the phase diagram of the 2d system using 
the calculated phase boundaries, eq. (18) and eq. (19), 
together with the data obtained from the molecular dy- 
namics simulations. Both the lattice-like and ring struc- 
ture can coexist in the area which is limited by these phase 
boundaries until the intersection is reached. If the temper- 
ature exceeds this boundary, only one state can survive de- 
pending on the field, leading to the loss of hysteresis. The 
pure ring phase is limited by the phase boundaries until T c 
is reached. If the critical temperature T c is reached at zero 
field, the rings dissolve and the particles enter the param- 
agnetic phase. For fields larger than (T) , the particles 
are unbonded, whereas the lattice-like structure is sustain- 
able only at low temperatures, T -C T c . The lattice-like 
structure gets lost due to the growing kinetic energy of the 
particles with increasing temperature. 

The phase diagram of the 3d system is depicted in fig. 6. 
We find the pure chain phase and a mixed phase where 
both rings and chains can coexist. The pure ring phase is, 
in contrast to the 2d system, limited to B = 0. Therefore, 
the area in which we can observe hysteresis is larger than 
in 2d systems. At high temperatures, T > T c , the system 
enters the paramagnetic phase, where the particles are 
unbonded. However, if the field is switched on the chain 
phase is again stabilized, because all particles are forced to 
point into the same direction. This makes an arrangement 
parallel to the field more favorable. At low temperatures, 
T < T c , and high fields, B > Bj, the particles can arrange 
only in chains. If the field drops to zero, the chain closes 



to rings provided that this bending is affordable. Other- 
wise, a short chain will turn into the direction of the field 
without closing to a ring. In this case, the system can go 
from the upper to the lower mixed phase without entering 
the pure ring phase. This behavior depends on the length 
of the chain because it is more favorable to bend longer 
chains to closed rings than shorter ones. This shifts the 
starting point of the pure ring phase from T = to higher 
temperatures. 

In this paper we have derived a phase diagram of dipolar 
nanoparticles as a function of particle size and tempera- 
ture, explaining the transition from magnetic chains and 
rings to a magnetic gas. Furthermore, we have calculated 
the phase diagrams for two and three dimensional systems 
as a function of temperature and external magnetic fields 
by molecular dynamics simulations. The different regions 
of stability explain qualitatively the collective behavior of 
magnetic nanoparticles which has been observed in several 
experiments. 
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